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We report measurements of microwave transmission over the first five Mie resonances of alumina 
spheres randomly positioned in a waveguide. Though precipitous drops in transmission and sharp 
peaks in the photon transit time are found near all resonances, measurements of transmission fluc- 
tuations show that localization occurs only in a narrow frequency window above the first resonance. 
There the drop in the photon density of states is found to be more pronounced than the fall in the 
photon transit time, leading to a minimum in the Thouless number. 



PACS numbers: 42.70.Qs, 42.25.Dd, 46.65. +g 

Photon localization has been of intense interest since 
it was first proposed as the analog of electron localiza- 
tion . Accurate measurements of the statistics of the 
electromagnetic field and intensity and of the total trans- 
mission in random samples have deepened our under- 
standing of mesoscopic physics and localization. How- 
ever, the very possibility of photon localization in three- 
dimensional systems without substantial crystalline order 
has been called into question by the barriers to achiev- 
ing and detecting electromagnetic localization. Unlike 
electrons that can be trapped by the Coulomb interac- 
tion at atomic sites, photons are not bound by individ- 
ual particles. They are not strongly scattered by par- 
ticles either, except at Mie resonances where the scat- 
tering cross-section can considerably exceed the geomet- 
ric cross-section. To preserve the strength of scattering, 
however, particles must be separated by at least the di- 
ameter of the scattering cross-section, requiring parti- 
cle separation that are larger than the wavelength A. It 
would then appear that the loffe-Regel-like condition for 
localization in three dimensions Q, isc < A/27r, where Isc 
is the scattering mean free path, cannot be satisfied by 
collections of independent resonant scatterers Mea- 
surements of exponential scaling of transmission j|, ||, ||] 
have not definitively established localization because the 
possible presence of absorption could also lead to expo- 
nential decay of intensity within a sample. However, re- 
cent measurements of coherent backscattering in macrop- 
orous GaP networks |^] along with theoretical predictions 
1^ suggest that the approach to localization can be ob- 
served in the rounded backscattering peak from weakly 
absorbing samples. Also recently, the variance of rela- 
tive fluctuations has been shown to provide a decisive 
test for localization, even in the presence of strong ab- 
sorption This provides a sure guide in the search for 
localization and raises the possibility that measurements 
of key scattering parameters can be used to sort out the 
precise material and structural characteristics that may 
edge samples towards and potentially across the localiza- 
tion threshold. 

In this Letter, we report microwave measurements of 
the frequency variation of three key localization param- 



eters in quasi-lD samples composed of randomly posi- 
tioned alumina spheres in an open copper tube. These 
arc the Thouless number |l^, which is the ratio of the 
width to the spacing between quasi-states of a random 
medium, 5 = (5z//Ai/, the dimensionless conductance g 
|]lT| , and a parameter g' Q , which represents the inverse 
of the variance of transmission normalized to its ensemble 
average value. These parameters capture, respectively, 
the relation to localization of average dynamics, average 
static transmission, and static fluctuations. Their mea- 
surement makes it possible to identify a narrow window 
of localization above the flrst Mie resonance and to probe 
a constellation of factors that foster localization, includ- 
ing size, concentration, and structural correlation. 

We report the first direct measurement of the Thou- 
less number. The Thouless criterion identifies the onset 
of localization as the point at which the level width Su 
becomes smaller than the level spacing Aj/ Beyond 
this point, transport is inhibited since modes in different 
blocks of a sample do not overlap. The level width 5v is 
identified with the field correlation frequency ||l^ . This is 
inversely related to the width of the time-of-flight distri- 
bution of the transmitted wave and to the average transit 
time through the medium, r p^ . The level spacing Av 
is the inverse of the density of states, Av = 1/N(y), thus 
giving 5 = 5i'N{v). The hurdles to localizing radiation, 
even in resonant media, can be seen from a dynamical 
perspective. The precipitous drop in 8v expected at reso- 
nance, reflecting the sharp peak in t, does not necessarily 
result in a diminished value of 5 because it is countered 
by an increasing value of Niy). 

The dimensionless conductance g is given by the sum, 
g = T^abTab 10 , where Tab is the transmission coefficient 
between incident mode a and output mode 6. In elec- 
tronic systems, g can be obtained from the measurement 
of the conductance, G = {e'^/h)g, whereas for classical 
waves it can be found from the measurement of the total 
transmission Ta — TibTab, since g — N{Ta), where N is 
the number of transverse modes given by = Ak'^ /2tt, 
where A is the area of the cross-section of a sample and 
k is the wave number, and the angle brackets denote the 
ensemble average. For diffusive waves in nonabsorbing 
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samples of length L, g = Ni/L, where i is the trans- 
port mean free path. The parameter g is essentially the 
number of diffusing channels in a sample and so the lo- 
calization threshold is reached when its value falls below 
unity [^. This is consistent with the Thouless criterion 
for localization, since for diffusive waves, g — 6 > I pl] |. 

In the diffusive limit and in the absence of absorp- 
tion, the variance of the total transmission normalized to 
its ensemble average is given by var{sa) — '^./'ig, where 
Sa = Ta/{Ta) 0, |l|l. Tliis allows us to construct a 
measure of fluctuations, g' = 2/3var{sa), which reduces 
to g in the above limit. The localization transition should 
occur when 5' ~ 1 [1). Since var{sa) and var{sab), where 
Sab = Tab/ (Tab), are related by var{sa) = ^{var{sab) - 1) 
0, we are able to evaluate g' either from measure- 
ments of total transmission, using an integrating sphere 
psf , or from local measurements of intensity ||] . 

Measurements of field transmission spectra are car- 
ried out in low-density collections of 0.95-cm-diameter 
alumina spheres (99.95% AI2O3) with use of a Hewlett- 
Packard 8772C network vector analyzer. Low densi- 
ties are produced by embedding the alumina spheres 
in Styrofoam spheres of low refractive index. The in- 
dex of refraction n of the alumina spheres is inferred 
from a comparison of the measured extinction coeffi- 
cient of the coherent intensity Ic in 3-D samples with 
an alumina volume fraction / = 0.010 with the cal- 
culated scattering cross-section asc for isolated alumina 
spheres. In the independent-scattering approximation. 
In 7c = —nsc<7scL, where is the number density of the 
alumina scatterers. The natural logarithm of Ic is com- 
pared to calculations of asc with n = 3.14 in Fig. 1. This 
gives the closest correspondence between the positions of 
the Mie resonances and those of the dips in In Ic . 

Measurements of microwave propagation and localiza- 
tion are made in quasi- ID samples with an alumina vol- 
ume fraction of / = 0.068 contained in a 7.3-cm-diameter 
copper tube of L = 80 cm and L = 12 cm. The alumina 
spheres are put at centers of 1.9-cm-diameter Styrofoam 
spheres. An ensemble of 5,000 statistically equivalent 
random samples is created by rotating the tube between 
successive field measurements. The incident mode a is 
parallel to the sample axis. The outgoing field is detected 
at a point on the output surface of the sample. For each 
sample configuration, the field spectrum yields the corre- 
sponding frequency variation of the transmitted intensity 
lab- The ensemble average (lab) for L = 80 cm samples, 
shown in Fig. 2a, exhibits distinct dips near each of the 
Mie resonances. The resonant character of scattering is 
further indicated by the sharp peaks in the spectrum of 
the average photon dwell time, seen in Fig. 2b. The tran- 
sit time is given by t = {sabdf/fab/duj), where (pab is the 
phase of the field and lu is the angular frequency [ p^ . 
Since low transmission can be due to absorption and 
long dwell time can be associated with microstructure 
resonances, these indicators do not provide a definitive 



measure of the closeness to the localization threshold. 
This can be obtained, however, from the measurement 
of var{sab), shown in Fig. 2c. For diffusive waves obey- 
ing Rayleigh statistics, var{sab) = 1- Lower values below 
8.5 GHz indicate a significant ballistic component in the 
transmitted field, while higher values indicate the pres- 
ence of substantial long-range correlation. The horizontal 
dashed line in Fig. 2c represents the localization thresh- 
old, var{sab) = 7/3, which corresponds to var{sa) = 2/3 
or, equivalently, g' — 1. The localization threshold in this 
quasi- ID sample is crossed in a narrow frequency range 
above the first Mie resonance. Large fluctuations in mea- 
surements of var(sab) above the localization threshold, 
seen in Fig. 2c, are a result of the greatly enhanced tail 
of the intensity distribution for localized radiation 

In order to take a closer look at the localization region, 
the data in Fig. 2 are superimposed and displayed using 
an expanded scale in Fig. 3a. The shift of extrema in 
(lab) and in var(sab) to frequencies above the maximum 
in T indicates that localization occurs above the resonant 
peak. Since localization occurs within a trough of (5 in a 
frequency range in which t is smaller than its peak value 
by a factor of three, N{i^) must drop from its peak value 
on resonance by an even greater factor. We also expect 
that in the absence of absorption, g' ~ g, and that the 
peak in var(sab) should occur at the same frequency as 
the dip in (lab) rather than at higher frequency as seen in 
Fig. 3a. The shift observed in the corresponding spectra 
in Fig. 3a is principally a result of the peak in photon 
dwell time, which increases absorption and suppresses 
(lab) more than var(sab) at the resonance. Enhanced 
absorption on resonance is also the cause of the drops in 
transmission at higher Mie resonances in Fig. 2a. 
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FIG. 1: Comparison of the scattering cross-section aac of a 
0.95-cm-diameter alumina sphere (solid) and the extinction 
of the coherent intensity Ic in a 3-D alumina sample with 
/ = 0.010 (circles). In calculations, the alumina sphere is 
assumed to be lossless and the refractive index is adjusted so 
that the positions of the Mie resonances are located at the dips 
in the experimental data. The curve shown is for n = 3.14. 
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To investigate the factors that inchne a sample towards 
locaUzation, we carry out the measurements presented 
in Fig. 3a in a shorter sample of L = 12 cm (Fig. 3b), 
so that £ < L < ^, I/a, where ^ = Ni is the localiza- 
tion length and La is the exponential absorption length. 
In this weakly absorbing sample, we find that transport 
is diffusive above 9.6 GHz, where 1 < var{sab) < 7/3 
(Fig. 3b), so that g, g', and S can be associated unam- 
biguously with the proximity to the localization thresh- 
old. The quantities g and g' are obtained from mea- 
surements of the total transmission T^, via expressions 
g — N{Ta) and g' = 2/3var{sa)- Good agreement be- 
tween g and g' is found, as seen in Fig. 4. Also notice 
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FIG. 2: (a) Average transmitted intensity {lab), (b) average 
photon dwell time r, and (c) var{sab) in a quasi-lD alumina 
sample of L = 80 cm and / = 0.068. The dashed line is the 
localization threshold. 
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FIG. 3: (color) (lat), t and var{sab) superimposed over the 
frequency range of the first Mie resonance: (a) L = 80 cm, 
(b) L = 12 cm. The dashed line is the localization threshold. 
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FIG. 4: (color) Localization parameters 5, g and g' over the 
frequency range of the first Mie resonance in an alumina sam- 
ple of L — 12 cm. The dashed line is the localization thresh- 
old. 



that, since the minimum oi g occurs above the resonance, 
the minimum of £ must also lie above that of £sc- 

To examine more closely the balance of competing ef- 
fects near the resonance, we measure Sv, Au, and d. The 
level width is determined from measurements of the 
field correlation function with frequency shift. We find 
that Si' shown in Fig. 5b is related to the average dif- 
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FIG. 5: (a) Part of a spectrum of d(j>ab/duj for a single con- 
figuration with tlie peaks marking the mode frequencies and 
(b) the level width 5v and the average level spacing Az/ in a 
L = 12 cm alumina sample. At all frequencies, 5v > Af. The 
dashed line is the inverse of N(v) = (^■kv'^ /vn) (vr, /ve ) V 
where Vp and ve are the phase M and the energy pcJ velocity, 
respectively, vo is speed of light in the effective medium, and 
V is the sample volume. 

fusion time (Fig. 3b), r = L'^/6D = l/Q5v. The level 
spacing Az^ is found by measuring transmission between 
two probes in the tube with copper end caps. In this 
case, line widths are narrowed and mode spacing can be 
measured (Fig. 5a). When determining Az/ both positive 
and negative peaks in spectra of d(j)ab/duj are accounted. 
The negative peaks appear because of interference be- 
tween two closely spaced modes. A comparison of Aj/ 
calculated for an ensemble of 360 sample configurations 
and Sv is shown in Fig. 5b. We find that, whereas dv and 
Av attain their minimum values at the resonance, their 
ratio is a minimum above the resonance. In the trough 
of S, where the measurements oi Av are most reliable 
by virtue of the lowest values of N{h'), (5 is in a good 
agreement with g and g' (Fig. 4). 

The increase above resonance in Ai/, obtained from 
a calculation of iV(z^) in the independent-scattering ap- 
proximation 1^, is well below the measured increase 
(Fig. 5b). This suggests that increasing the density of res- 
onant scatterers beyond the point at which they begin to 
interact suppresses the density of states and thereby fos- 
ters localization. A decrease in N{iy), on the other hand, 
was predicted by John pl| ] for nearly periodic samples. 
Indeed, for periodic samples exhibiting a complete pho- 
tonic band gap, the vanishing of iV(j/) automatically re- 
sults in localization within a pseudogap, when disorder is 
introduced. We find that, when the alumina spheres are 
displaced by half of their radius from the centers of the 
embedding Styrofoam spheres, var{sab) is only slightly 
suppressed. For a L = 55 cm alumina sample, var{sab) 



drops from 6.0 to 4.5 at its maximum. This small reduc- 
tion, for a considerable broadening of the radial corre- 
lation function, suggests that collective scattering from 
randomly positioned alumina spheres itself and not only 
their residual positional order is the chief reason for the 
sharp drop in N{i') above the first Mie resonance. 

In conclusion, we have measured the dimensionless 
conductance, the variance of transmission fluctuations, 
and the Thouless number in low-density collections of 
alumina spheres in a quasi-lD geometry. We find these 
are equivalent measures of localization for L < La,£,. 
For samples longer than La, only var(sa) serves as a 
reliable guide to localization, which is found in a nar- 
row frequency range above the first Mie resonance in the 
longer sample. Localization does not occur on resonance 
because the narrowing of the level width is offset by the 
enhancement of the density of states. However, a sharp 
drop in N^iy) above resonance leads to a minimum in 
S. The drop in N(h') and the consequent appearance of 
localization are the result of collective scattering. 
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